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Final  Report 

Award  Number  “N000 14-96- 1-0264” 


LONG-TERM  GOALS,  OBJECTIVES  AND  APPROACH 

The  long-term  goal  over  this  three  year  project  has  been  to  develop  computer  software 
needed  to  optimize  initial  conditions,  internal  parameters  and  external  parameters  for  the 
Harvard  primitive  equation  (PE)  model  in  order  to  produce  the  best  forecasts  in  an  arbitrary 
region.  This  new  tool  invokes  an  inverse  technique  to  fuse  all  available  data  types,  gathered 
non-synoptically,  with  optimized  model  dynamics.  The  technique  is  distinct  from  (and 
complementary  to)  the  optimal  interpolation  and  Kalman  filter  assimilation  strategies  now 
being  developed  and  used  at  Harvard  (e.g.,  Lermusiaux,  1997). 

The  scientific  objectives  of  this  research  include  answering  the  following  questions.  Can 
forecast  skill  in  a  highly  unstable  region  like  the  Iceland-Faeroe  Front  be  extended  to  7  days? 
Can  a  diagnostic  simulation  over  a  10-day  interval  in  that  region  include  all  the  data  in  an 
inverse  calculation,  or  is  it  too  nonlinear?  What  are  the  relative  impacts  of  the  various  data 
types  (CTD/XBT/XCTD  casts,  current  meters,  surface  drifters)  on  making  forecasts  in  this 
region?  The  technical  objectives  encompass  the  details  of  the  model  fitting  process.  How 
nonlinear  is  the  fit?  Can  the  nonlinearity  be  reduced  by  optimizing  large-scale  structure 
first?  How  much  data  can  be  fit  at  one  time?  Is  the  distribution  of  the  data  sufficient  to 
initialize  the  model?  Are  the  open  boundaries  causing  instabilities  in  the  model? 

Real-time  ocean  forecasting  involves  assembling  an  initial  state  which  often  requires  merging 
many  data  types  that  are  usually  gathered  over  non-synoptic  intervals.  Furthermore, 
dynamical  ocean  forecast  models  still  require  improvements  in  their  physics  (including 
parameterizations).  We  are  addressing  these  two  issues  simultanously  in  applying  a 
standard  inverse  technique  to  the  Harvard  PE  ocean  model  in  the  context  of  an  unique 
dataset  in  the  Iceland-Faeroe  frontal  region.  The  computational  strategy  that  we  are 
developing  will  be  applicable  to  whatever  new  regions,  updated  PE  model  versions,  or 
varieties,  amounts  or  qualities  of  data  become  available.  The  Harvard  ocean  modeling  group 
provided  expertise  in  setting  up  and  modifying  the  model  for  our  purposes. 
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WORK  COMPLETED 


During  the  three-year  project,  we  prepared  the  data,  designed  the  model  domain,  developed 
the  iterative  process  for  minimizing  model-data  mismatch  and  applied  the  method  to  both  an 
identical-twin  experiment  and  the  data  from  the  1993  IFF  surveys.  The  results  of  our 
research  were  submitted  to  the  special  issue  of  Dyn.  Atm.  Oceans  dedicated  to  Prof.  A.  R. 
Robinson  (Miller  and  Comuelle,  1998). 

Miller,  A.  J.,  Lermusiaux,  P.  F.  G.  and  Poulain,  P.-M.,  1996.  A  topographic-Rossby  mode 
resonance  over  the  Iceland-Faeroe  Ridge.  J.  Phys.  Oceanogr.,  26:2735-2747. 

Miller,  A.  J.  and  B.  D.  Comuelle,  1998:  Forecasts  from  fits  of  frontal  fluctuations.  Dyn. 
Atmos.  Oceans,  sub  judice. 

RESULTS 

The  fitting  procedure  was  tested  and  applied  with  the  Harvard  PE  model  to  non-synoptic 
hydrographic  surveys  of  unstable  current  meandering  of  the  Iceland-Faeroe  Front.  The 
initial  conditions  (including  the  fields  outside  the  data  domain)  were  adjusted  with  large- 
scale  structure  functions  to  fit  the  observations,  with  no  additional  forcing  or  adjustment  of 
boundary  conditions  during  the  model  runs.  The  technique  was  first  tested  with  an  'identical 
twin’  predictability  experiments  that  showed  the  iterative  technique  can  fit  the  non-synoptic 
Initialization  Survey  data,  partly  correct  the  initial  condition  errors  and  allow  the  model  to 
move  closer  to  the  true  evolution.  However,  the  limited  non-synoptic  verification  data  was 
inadequate  for  providing  convincing  evidence  of  forecast  skill  even  for  short  2-day  forecasts 
which  are  known  to  have  skill  in  the  identical  twin  framework. 

The  PE  model  was  then  fit  to  the  observed  hydrographic  data  from  August  1993  in  several 
scenarios.  With  only  a  few  iterations,  typical  temperature  and  salinity  along-track  error 
variance  reductions  of  70-80%  were  achieved  relative  to  initializing  the  model  from  a  time- 
independent  objective  analysis.  Although  hindcast  skill  increased  by  applying  the  method, 
significant  quantitative  forecast  skill  was  not  achieved  probably  due  to  inadequate 
initialization  and  verification  data  (as  found  in  the  twin  experiments).  Qualitative  skill 
assessment  proved  necessary  to  distinguish  the  integrity  of  the  hindcasts  and  forecasts, 
especially  with  respect  to  the  occurrence  of  a  hammerhead  baroclinic  instability  of  the  IFF 
(Miller  et  al.,  1995;  Robinson  et  al.,  1996).  Since  the  model  was  successfully  fit  to  the 
hammerhead  instability,  the  incorrect  hammerhead  forecasts  from  fits  must  be  a 
consequence  of  (a)  inadequate  initialization  data,  (b)  high  sensitivity  to  initial  conditions,  (c) 
incorrect  model  physics  or  (d)  insufficient  structure  functions.  The  quantitative  forecast 
skill  shown  by  Robinson  et  al.  (1996),  who  used  an  optimal  interpolation  technique  in  real 
time,  is  comparable  to  what  was  obtained  here  for  the  hammerhead  when  the  model  was 
initialized  from  or  fit  to  the  antecedent  (Zig-Zag)  survey. 


The  main  difference  between  this  method  and  the  representer  method  of  Bennett  (1992)  is 
our  use  of  the  restricted  set  of  basis  functions.  An  adjoint  model  can  give  the  complete 
structure  of  the  data  sensitivity  in  the  forward  problem,  but  using  similar  smoothing 
assumptions  (i.e.  a  covariance  matrix  for  the  initialization  changes)  should  lead  to  similar 
solutions.  The  restriction  to  large-scale  perturbations  is  useful  for  linearized  methods  since 
the  small-scale  structure  of  the  sensitivity  should  be  less  reliable  than  the  large-scale 
structure,  due  to  nonlinearity.  For  example,  in  the  identical  twin  experiment  where  the  error 
was  spanned  by  50  structure  functions  the  linearized  iteration  did  not  produce  the  exact 
initial  condition.  A  similar  problem  would  arise  for  the  linearized  derivatives  made  by  an 
adjoint  model. 

The  model  run  from  the  optimized  initialization  constitutes  a  possible  dynamically 
consistent  scenario  explaining  some  of  the  variability  seen  in  the  IFF  observations  as  a 
meandering  of  the  front.  To  the  extent  that  the  model  is  accurate,  the  observations  have  been 
reconstructed  into  a  four-dimensional  picture  of  the  flow  field  in  the  area.  The  results  here 
set  the  stage  for  diagnostic  analyses  of  the  frontal  baroclinic  instabilities,  e.g.,  using  a 
primitive  equation  analogue  (under  development  by  Miller,  Lermusiaux,  Pinardi  and 
Robinson)  of  the  energy  and  vorticity  analysis  (EVA)  system  developed  by  Pinardi  and 
Robinson  (1986)  for  quasigeostrophic  dynamics  and  applied  by  Miller  et  al.  (1995)  to 
quasigeostrophically  modeled  IFF  variations.  The  development  of  the  EVA  system  was  also 
supported  by  ONR  under  a  supplement  to  the  main  proposal.  We  set  up  subroutines  for 
the  mapping  from  sigma  coordinates  to  z-levels  as  well  as  the  basic  structure  of  the  energy 
balances  (dervied  from  the  original  codes  of  Spall).  We  will  continue  to  develop  these 
routines  in  the  coming  years  and  test  them  in  California  Current  modeling  studies. 

A  correct  ocean  model  is  a  valuable  tool  for  the  interpolation  and  interpretation  of  data,  as 
well  as  practical  prediction  problems,  but  it  is  dependent  on  obtaining  adequate  data  for 
judging  model  quality.  The  techniques  tested  here  could  easily  accommodate  other  data 
types,  but  the  best  type  of  data  for  constraining  the  model  initializations  and  for  testing 
forecast  skill  would  appear  to  be  synoptic  hydrographic  surveys  as  could  be  obtained  from 
aerial  XBT  surveys.  For  this  August  1993  IFF  dataset  there  are  drifter  observations, 
current  meter  observations  and  a  satellite  SST  image  which  can  be  used  in  future 
applications  of  this  technique  to  improve  the  fits  and  forecasts.Giving  higher  vertical 
resolution  to  the  model  is  of  high  priority  in  improving  the  details  of  the  model  fit, 
especially  the  baroclinic  structure  of  the  hammerhead  instability.. 
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